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Abstract 

We have observed a negative differential conductance with singular gate and source-drain bias de- 
pendences in a phosphorus-doped silicon quantum dot. Its origin is discussed within the framework 
of weak localization. By measuring the current-voltage characteristics at different temperatures as 
well as simulating the tunneling rates dependences on energy, we demonstrate that the presence 
of shallow energy defects together with an enhancement of localization satisfactory explain our 
observations. Effects observed in magnetic fields are also discussed. 
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I. INTRODUCTION 



The consequences of localization and disorder on the transport properties, in particular, 
on the mobility of devices have been investigated extensively in bulk and in two-dimensional 
semiconductors [l, 2|. One of the most noticeable and known effect is hopping conduction 
where electrons are allowed to hop between impurities sites, leading to a decrease in the 
material conductivity with temperature. This was first theoretically described by Abrahams 
[sl to explain the transport properties of random media. In this case, the boundary between 
ballistic and diffusive transport is well represented spatially by an ellipsoid which major and 
minor axis are given respectively by the hopping and localization lengths. Unless screened 
by a metal gate or short-range disorder, electron-electron interaction is often competing 
with localization, leading to a modification of the conductivity at low temperature ^ and 
eventually to a metal to insulator transition jsj]. In nanostructures and, in particular in 
quantum dots, the localized region is defined by the device geometry, principally the number, 
the shape, the nature and the orientation of interfaces. Due to dielectric screening, the 
effective Bohr radius at trapping sites at the edges of the structure can significantly increase 
and electron localization becomes possible. Such phenomenon have been already observed at 
the periphery of dots fabricated within an AlGaAs-GaAs heterostructure {g] and Coulomb 
blockade has been found to be affected via a change in the dot capacitance. 

In doped devices, like phosphorus-doped silicon single electron transistors, the localized 
region has a more profound effect owing to the presence of dopants at the edge of the 
structure and so, in close proximity to the Si-Si02 interface j?!. In presence of Coulomb 
interaction, charge rearrangement is possible and becomes an effective mean of reducing the 
charging energy. It may also be responsible for the observed reduction in the effective dot 
dimension due to electron trapping at the edge of the structure and the subsequent change 
in the electrostatic potential. 

Despite this, the outer region of the dot is not the only place where localization has 
a non-negligible influence. Randomness in dopant distribution may lead to traps being 
present at or close to the tunnel barrier. As a result, the transmission probability, and so, 
the device conductivity, is expected to be significantly affected. The situation where a single 
dopant is unintentionally located at the barrier has already been investigated js, 9| in silicon 
nanostructures. However, most studies report on resonant features following the alignment 



of the quantum dot levels with the Dq or D_ states of the impurity 10| , a system that drives 



interest as it could be applied to spin filtering 



11| in spin qubit architectures [1 



In this article, we have used a highly doped silicon single electron transistor (SET) and 
show that an additional localization process may lie at the tunnel barrier. Caused by 
natural defects or more generally, by an asymmetry in the barrier shape, the effect can 
enhance significantly the electron trapping mechanism. In the first section, we discuss 
the observation of a negative differential conductance as well as its dependence on gate 
voltage. We then describe the localization process through the temperature dependence of 
the conductivity and simulations. The analysis of experimental results allows distinguishing 
between static and dynamic effects in the tunnel barriers. In a second paragraph, we describe 
and interpret the magnetic field and temperature dependences in the low temperature regime 
before conclusions are drawn on the origin and potential implication of this effect. 



II. DEVICE AND EXPERIMENTAL SETUP 



Devices were fabricated from a silicon-on-insulator (SOI) substrate with a 40 nm-thick 
silicon layer that was subsequently implanted with phosphorus to obtain a dopant density 
of 3x10^^ cm~^ after thermal anneal. At such a doping density, the largest areas of the 
contact leads remain metallic whereas high density defect surfaces partially deplete the 
nanostructures and locally increase the localization at its edges The quantum dot and 
its controlled in-plane gate were defined by electron beam lithography during a single 
step process. The dot diameter was then reduced to about 120 nm after thermal oxidation, 
leading to the formation of a 20 nm oxide layer (Inset Fig. 1). All measurements were taken 
in an Oxford Instrument Heliox"^^ cryostat that was fitted with a 7 T magnet and whose 
base temperature was 300 mK. The current /so was measured using a Hewlett Packard 
3458A multimeter with a 10^ V/A amplifier on battery. The temperature of the device was 
controlled by an Oxford Instrument Intelligent Temperature Controller (ITC4) from 300 mK 
up to 150 K using Ru02 and Cernox calibrated thermometers. 
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FIG. 1. Position of the NDC in VgD and at 300 mK. The inset shows a Scanning Electron 
Microscope (SEM) image of the device with source (S), drain (D) contacts as well as the SET gate 
Vg. The bar scale represents 200 nm. 



III. EXPERIMENTAL RESULTS 



A. Negative differential conductance 



The device is characterized by regularly spaced Coulomb diamonds and so, it operates 
as a single quantum dot at low temperature. This will not be discussed here, but one may 



report to previous observations [15|. However, a negative differential conductance (NDC) 
can clearly be noticed in the region of negative source-drain bias V^-d (Fig. 1). Such feature 
could be found in single quantum dot where additional parallel channels via edge states 



are opened for tunneling or when trap-assisted tunneling is present 13| or in the case of 
spin dependent transport [l^. Here, the position of the NDC clearly follows the shape of 
the Coulomb diamonds, at a fixed distance from their edges. This rather suggests that the 
existence of the NDC is associated with tunneling events in the quantum dot. Previous 
studies on the same quantum dot have shown that such a NDC was offsetting the electron 



tunneling detection point of a nearby SET 15[ and that an asymmetry in the tunneling rates 



could explain the observed phenomenon. Still, this observation alone does not enlighten on 
the fundamental electronic process responsible for such asymmetry in the tunneling rates 
and so, additional experiments are required. 
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FIG. 2. Conductivity in the inelastic cotunneling regime at different values of Vso- 



B. Localization effects in the tunnel barrier 



1. Barrier asymmetry 



The asymmetry of the tunnel barrier rates can be assessed more directly by deter- 
mining the dependence of the multi-electron tunneling rates for each barrier 16[ /"s, d 
on VsD- Extracting separately the values for Fg and Fj^ from the total tunneling rate 
-^sd/c = FsF^/ {Fg + Td) and the dependence of the current on temperature is generally 
challenging. This is due to the presence of both elastic and inelastic processes and their 
complex dependences on source-drain and gate biases. 

However, it is possible to estimate the product of the tunneling rates by observing the 
quadratic dependence of the conductivity in temperature below 3 K, a behavior that is char- 
acteristic of inelastic cotunneling (Fig. 2). Indeed, within this regime the conductivity 
is given by 



2he' 
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(1) 



where Ec is the dot charging energy, AF = —eaoVsu — eOgVg the change in the free 
energy of the system on the tunneling of an electron in or out the quantum dot, and a^, 
the level arms of the drain and gate respectively. 



On the other hand, it is possible to estimate the tunnehng rates by taking into account the 
tunnel probabilities of an electron through a barrier at a specific energy level En = NEq + Eq 
in the quantum dot. The use of the parabolic barrier approximation 1^ allows extending 
the range of validity of the tunnel rates to large V^d. Following, a method similar to the one 
used for square barriers 19|, we found : 



To oc ^^^exp [a {AF + En) + /?] /(eanV^so, Te 



(2) 



f{E,T) is the Fermi distribution of the drain reservoir, T^, is the electron temperature, 
AE the level spacing of the dot. Coefficients a, (3 are given respectively by : 



a = W 

2 [ 
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/3 = 

' 2 \h^v 



1/2 



(4) 



where m* = 0.19mc_ is the effective mass of silicon, V and W respectively the height and 
width of the tunnel barrier (defined at the Fermi energy of the reservoir). 

The value of FsF^ is then calculated by summing over all possible tunneling events, 
including elastic and inelastic contributions for the different energy levels of the quantum 
dot (ground and first excited states). The behavior of /"s-^d is well reproduced for both 
positive and negative Vgn (especially at small biases) (Fig. 3) if one takes into account the 
dependence of the electron temperature on source-drain bias due to acoustic phonons 2ot : 
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where Tq is the phonon bath temperature and the normalization constant Vq = 1 mV. 

The drain and gate level arms can be extracted from Coulomb diamonds and we find 
ttg = 0.048 and ttgn = 0.44. From best fits, we obtain Ws = 8.0 nm and Wu = 13.3 nm 
for the source and drain barrier widths and, = 5.4 meV and Vq = 11.8 meV respectively 
for the corresponding barrier heights. The values for the barrier widths are compatible with 
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FIG. 3. Variation of the product of the tunnehng rate of the source {Fs) and the drain (-Td) with 
source-drain bias and corresponding simulations for parabohc barriers (dotted hues) and in the 
case of energy dependent barrier width (sohd hne). 
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FIG. 4. Current-voltage characteristics before (solid line) and after thermal cycle to 150 K (dotted 
line) and to room temperature (dashed line). 



the ones obtained from SEM observations and suggest the existence of an asymmetry in the 
device structure, and so, an asymmetry in the tunneling rates. It is interesting to notice 
that the potential barriers are small implying parallel conductive path at relatively small 
bias values, an effect explaining the additional conductivity background for | V^^ |> lOmV. 

There is, however, a noticeable deviation for Vso < — 10 mV between the simulated curve 



using the parabolic approximation for the potential barrier and the experimental data (Fig. 
3). This indicates that the tunnel rate is substantially decreased below this value. Because 
the tunneling time in a single square barrier is given by Tg ~ sinh (201^3), where a is the 
attenuation constant, the potential barrier has to be substantially wider at energies lower 
than 5.3 meV (below the Fermi energy of the source lead when Vsd=0). This effect is well 
confirmed by simulations if electron tunneling is significantly decreased for energy levels 
lying below -AEc- To investigate the origin of such barrier shape, and in particular, to 
distinguish between static (geometry) and dynamic (localization) effects, we proceeded with 
a series of thermal cycles (Fig. 4). 

After raising the device temperature to 150 K and slowly cooling it back to 300 mK, no 
NDC is observed. Instead, a plateau is present in the current- volt age characteristics. Follow- 
ing a different cooldown from room temperature, no NDC nor plateau could be observed. 
These results point towards the presence of trapped electrons that have been thermally 
activated. Further thermal cycles up to room temperature did not change the shape of 
the current voltage characteristics and only an inflexion remains (arrow in Fig. 4). This 
indicates that both a geometrical asymmetry in the barrier and electron trapping during 
cooldown are likely to be responsible for the NDC. 

2. Enhancement of localization 

Additional information on the trapping mechanism is obtained by analyzing the tem- 
perature dependence of the conductivity in the high temperature regime. To this end, we 
proceeded to a thermal cycle between 300 mK and 150 K and measured the variation in the 
height of a Coulomb peak in temperature (Fig. 5). The behavior is significantly different 
depending on whether the temperature is increased (squares) or decreased (circles). In the 
first case, the conductivity is limited by cotunneling at low temperature, as shown before 
(Fig. 2), then follows an activation behavior between 5 and 15 K, with an associated energy 
El = 1.39 meV. The most noticeable feature is the abrupt decrease in conductivity at To ~ 
15 K and the presence of a new activation process with an energy value E2 = 3.5 meV above 
this temperature. In the second case, i. e. when the temperature is gradually decreased, 
the same activation process E2 is present in the high temperature range but the process 
El has disappeared. This suggests that E2 is an intrinsic characteristic of the device, and 
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FIG. 5. fj (r) with activated (solid lines) and cotunneling (CT) regimes when increasing and 
decreasing temperature from 3 to 150 K as shown by arrows. The dotted line shows the region of 
trap ionization associated with the electron delocalization temperature Tq. Vso = 1-5 mV and 
= -2.27 V. 

so, corresponds to an activation of donors in silicon, whereas Ei is associated with the way 
electrons are localized to impurity states during the cooldown. Such an irreversible and 
abrupt delocalization of traps is supported by the fact that ksTo ~ Ei. For single donors in 
silicon or at low concentration, the ionization energy is expected to be close to 45 meV. This 
value, however, decreases substantially with the donor separation, and, at a concentration of 



2l|. As 



~ 10^^ cm~^ similar to the one used in these devices, it is presumably close to zero 
a consequence, E2 is likely to be associated with phosphorus donors at the Si/Si02 interface 
where localization effects may be more important. This increase of localization arises from 
the sensitivity of the ionization energy to the dielectric environment. Indeed, for degenerate 
silicon, we have |22| : 



where z is the distance to the Si02-Si interface and Eq the dielectric constant in free space. 
Indeed for 2; ~ 9 nm, we obtain E2 ~ 3.5 meV. 
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FIG. 6. Temperature dependence of the conductivity in the low temperature regime. A change in 
behavior is clearly noticeable at around 1.73 K. Ec is the charging energy of the quantum dot. 

Finally, we notice that the NDC disappears irreversibly above 17 K, thus linking the 
electron localization in the barrier to the presence of an NDC. Still, this process alone 
cannot explain the observed phenomenon as a whole. Instead a combination of structural 
defects in the barrier and localization enhancement due to electron trapping may explain 
the behaviors observed in temperature and gate voltage dependences. Structural defects in 
the barrier may arise from dislocations or the presence of impurities, silicon dangling bond 
defects or the formation of a phosphorus clusters that localize electrons that may form either 
during the implantation or the annealing stage. 

After thermal cycling, the low temperature conductivity is higher, implying that the 
strength of localization has been decreased successfully. 



C. Low temperature and magnetic field dependences 

1. Temperature effects 

Despite the high doping concentration, the silicon quantum dot shows an insulating 
behavior down to 300 mK, for all gate voltages i.e. the conductivity decreases with 
temperature (Figs. 2 and 5). However, for a small range of source-drain biases Vgo close in 
value to the the NDC position and in the absence of magnetic field, we observe an increase 
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FIG. 7. (a) Variation of the position of the NDC in gate voltages with temperature, (b) Evolution 
of the shape and position of the NDC at 0.3 K (solid line), 1.5 K (dotted line) and 3 K (dashed 
line). The apparent shift towards positive Vsd is an effect of the displacement of the NDC in gate 
voltage, as described in the text. 

in conductivity a below = 1.73 K (~ 149 /leV) (Fig. 6). To understand the origin of 
the effect, we have measured the dependence of the conductivity on Vsd and V^, similarly 
to Figure 1, but at different temperatures up to 3 K . We find that the NDC position (Vg^ 
and V*) is significantly affected by the temperature (Fig. 7a). Also, the NDC evolves into 
a plateau due to the presence of additional tunneling mechanisms above 1.5 K, a value close 
to Tc (Fig. 7b). 

By raising the temperature of the device, the conductivity increases via cotunneling and 
the energy levels get broader, leading to a loss in the NDC visibility. The second effect of 
thermal broadening is a slight displacement of the NDC in source-drain bias of the order 
of 3.53 AjbT [23|. This corresponds to the thermal broadening of a Coulomb peak in the 
classical and high temperature regime, so the expected shift at 1.2 K is about 830 fiV, as 
measured (Fig. 7). The increase of temperature also modifies the tunneling rates via Eqs. 1 
and 2, so that the fractional number of electrons in the quantum dot n becomes temperature 
dependent and 

n = — ^ ~ 1 - 7exp {5E/k, T) (7) 

D + S 

where 7 is a constant depending on the device characteristics and 5E a function of Vst, 
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FIG. 8. a) NDC position at different magnetic fields (field number in red in inset) after thermal 
cycle, b) I^u (Vsd) at (solid line) and 5 T (dotted line) before thermal cycle and at a fixed Vg. 
The relative similar positions of the NDC is an effect of the dependence of the NDC position on 
B. c) Sketch of the effect of magnetic field on the density of states of the localized states located 
in the barrier. 

and Vg. 

Consequently, Coulomb diamonds shift along Vg with temperature, so is the position of 
the NDC, as observed experimentally. Since Figure 6 is obtained at fixed values of Vg and 
VgD, the variation da/dT reflects indeed the displacement of Vg-^ and V* in the map of Figure 
1, and consequently, the variation of the conductivity a (V^, Vsd) at low temperature. 

2. Magnetic field effects 

Similarly to temperature effects, the application of a magnetic field B perpendicular to the 
substrate modifies both the position of the NDC and its visibility (Fig. 8a). One may expect 
an increase in the localization strength at high 5, and so, a more pronounced NDC. However, 
only an inflexion in the current voltage characteristic is visible (Fig. 8b). Unlike thermal 
cycles, this effect is reversible. The effect of the field does reduce the effective diameter of 
the quantum dot by localizing more strongly the electrons at the Si-Si02 interface. It also 
affects the capacitances values. However, this is the density of states (DOS) near the source 
barrier which is predominantly affected. The number of localized states in the barrier is 
conserved but the spread in energy level is reduced (band tails in the DOS disappearing) 
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(Fig. 8c), so the NDC is smoothed. Such a modification of the DOS with B also explains that 
the tunneling rates are locally increased when tunneling takes place. The effective electron 
number in the dot is then reduced and the NDC position in is shifted accordingly. The 
combination between an increase in the charging energy and the disappearance of the DOS 
tails justifies the magnitude of the shift observed along Vgo. 



IV. CONCLUSIONS 



We have shown that the NDC that we observed could be explained by the presence of 
defects in one of the tunnel barrier and by an increase in the localization in that region due to 
electron trapping during the cooldown of the sample. Within this description, the change in 
the sign of da/dT results from the combination of the presence of an NDC region at negative 
source-drain bias consecutive to localization effects and the displacement of its position in 
Vg and Vsu with temperature, an effect which is partly due to thermal broadening and a 
modification of the tunneling rates. The application of a perpendicular magnetic field mainly 
affects the DOS, leading to similar effects to the temperature, including a displacement of 
the NDC with respect to gate and source-drain biases and the disappearance of the NDC at 
high field. These effects may have implications on the electrostatic control over the tunnel 
barriers as well as the spin-dependent electron tunneling if such a method is used for reading 
out the spin states of an electron. 
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